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! Abstract 

a^ 

The Wheeler-DeWitt equation is based on the use of canonical quantization 
rules that may be inconsistent for constrained dynamical systems, such as 
minisuperspaces subject to Einstein's equations. The resulting quantum dy- 
namics has no classical limit and it suffers from the infamous "problem of 
time." In this article, it is shown how a dynamical time (an internal "clock") 
can be constructed by means of a Hamilton- Jacobi formalism, and then used 



m 
(N 



' for a consistent canonical quantization, with the correct classical limit. 

O ■ 1 Introduction 
O 

I Classical field theories describe physical phenomena by means of field variables subject 

^ ' to partial differential equations. It often happens that the number of field variables exceeds 
^ • that of the physical degrees of freedom: there is no unambiguous way of prescribing 
I . the values of the field variables that correspond to a given physical situation. In the 
1^'. mathematical structure of the theory, this property is reflected by the existence of a 
gauge group, that allows transformations of the field variables while the physical situation 
'■ remains unchanged. 

' The peculiar feature of Einstein's theory of gravitation, which sets it quite apart from 

^ ordinary field theories, is that its gauge group consists of arbitrary distortions of the space- 

time coordinates, and thus cannot be disentangled from the structure of space-time itself. 
In particular, the time evolution of the gravitational field is locally indistinguishable from 
a gauge transformation — namely a local distortion of the space-time coordinates. As a 
consequence, the Hamiltonian density Ti, which generates the time evolution of the field 
variables, is weakly equal to zero [1]: namely, although 7i is a nontrivial function of 
the field variables (so that it can generate a nontrivial evolution), its numerical value is 
constrained to vanish. 

This Hamiltonian constraint does not cause any difficulty in the classical canonical 
theory. The numerical value of the Hamiltonian is only an initial value constraint. It is 
the functional form of the Hamiltonian that is needed for deriving the equations of motion. 
In quantum theory, however, if the gravitational field equations are quantized according 
to the standard canonical rule, namely vr'^"' = —ihS/6gmn, the resulting Wheeler-DeWitt 
equation [2, 3] leads to a dilemma known as "the problem of time." The difficulty is that, 
when the constraint = is imposed on the state vector ^, the latter is "frozen." There 
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cannot be wavepackets moving along classical trajectories, in accordance with Ehrenfest's 
theorem [4]. That is, the quantum equations do not lead to the expected classical limit. 

The problem was investigated long ago by Arnowitt, Deser, and Misner [5] who showed 
that there is an infinite number of possible coordinate conditions that may be used to 
put the theory in canonical form. The imposition of these coordinate conditions is equiv- 
alent to the introduction of "intrinsic" coordinates, defined by the dynamical variables 
of the physical system. The ADM method [5] and similar ones [6] provide, in principle, 
a completely general solution to the problem. Unfortunately, it is difficult to actually 
implement such a solution for a given, specific, physical situation. 

In particular, many authors have been interested in the properties of highly sym- 
metric cosmological models, for which there is a reasonable hope of obtaining an explicit 
solution. The trouble is that these symmetric situations have fewer (if any) dynamical 
degrees of freedom to which the ADM conditions can be applied, and the problem of time 
arises again. It would be impossible to mention here all the attempts that were made to 
solve that problem. Only a few randomly chosen references are listed below [7-11], with 
apologies to the authors of many similar works. 

Why is general relativity special? The reason is the totalitarian nature of Einstein's 
theory: all physical phenomena are coupled to gravitation, for the simple reason that all 
phenomena occur in space-time, and the properties of space-time are determined by the 
gravitational field. Any stresses, or any other physical forces, are themselves sources of 
the gravitational field, by virtue of the Einstein equations. (Electromagnetic theory, as a 
counterexample, is compatible with the existence of forces of non-electromagnetic nature, 
that have no electromagnetic field of their own, and can cause electric charges to move in 
arbitrary ways.) 

On the other hand, quantum theory, unlike general relativity, is not a "theory of every- 
thing" [12]. Its mathematical formalism can be given a consistent physical interpretation 
only by arbitrarily dividing the physical world into two parts: the system under study, 
represented by vectors and operators in a Hilbert space, and the observer (and the rest 
of the world), for which a classical description is used. This point was emphasized long 
ago by Bohr [13]: 

The necessity of discriminating in each experimental arrangement between 
those parts of the physical system which are to be treated as measuring in- 
struments and those which constitute the objects under investigation may be 
said to form a principal distinction between classical and quantum-mechanical 
description of physical phenomena . . . The place within each measuring pro- 
cedure where this discrimination is made is largely a matter of convenience. 

The consistency of this hybrid quantum-classical formalism can formally be proved, under 
suitably restrictive assumptions on the properties of the classical world [14]. That is, as 
foreseen by Bohr, the precise location of the boundary between the classical and quantum 
parts of the system is irrelevant for well posed problems. 

General relativity and quantum theory therefore appear to be fundamentally incompat- 
ible. Nevertheless, as will be shown in this article, they can be combined in a consistent 
way, by a careful choice of the dynamical variables. This is shown below by means of 
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two simple examples, with a few degrees of freedom. Each example starts by specifying 
a Lagrangian. Indeed it is known that canonical commutation relations are compatible 
with specified equations of motion only if the latter are equivalent to the Euler-Lagrange 
equations derived from some Lagrangian [15]. In both examples, the Lagrangians are 
chosen in such a way that the resulting dynamics are afflicted by the infamous "problem 
of time," just as in canonical quantum gravity. Yet, canonical quantization is possible, 
provided that one degree of freedom is kept classical, so that it can be used as a clock. 

The general method for solving this type of problem, for an arbitrary Lagrangian, is 
briefly discussed. It involves the solution of a flrst order partial differential equation, 
similar to the Hamilton- Jacobi equation, but with a very different physical meaning. 



2 A simple example of constrained dynamics 



As our flrst example, consider a dynamical system with three degrees of freedom, x, y, z, 
and with a Lagrangian 



2 I z 



zx 



z 



The canonical momenta are px 
thus are 



(1) 



xj z, py = —y/z, and pz = 0. The equations of motion 



Px = d{x/z)/dt = dL/dx = —zx, 
Py = d{-y/z)/dt = dL/dy = zy, 



Pz = = dL/dz = i {-pI -x^+pI + y") 



(2) 
(3) 
(4) 



Equations of motion with a similar behavior occur in a cosmological model with a homo- 
geneous scalar fleld, as in the Weinberg-Salam theory, but without a Higgs potential [16]. 

It is now convenient to introduce an auxiliary time, r, by means of dr = zdt, so 
that Px = dx/dr and Py = —dy/dr. Together with this new time, we also have a new 
Lagrangian, given by Lrdr = Ldt, so that the action remains the same. That is. 



Lt- = L/z 



dy^ 
dr 



The corresponding Hamiltonian is 

Hr = \{pl + ^'-pl-y')^ 



(5) 



(6) 



and it is easy to derive from it the equations of motion of the two harmonic oscillators, 
X and y. It follows from Eq. that Hr = 0. This causes no difficulty at the classical 
level. The Hamilton equations of motion are derived from the functional form of the 
Hamiltonian, irrespective of its numerical value, and they are equivalent to Eqs. (|^) and 
above. 
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Trouble arises, however, if we attempt to quantize such a system by introducing a wave 
function ip{x, y) that satisfies Hrip=0. Separation of variables readily leads to the general 
solution 



1p = ^2*^^ Un{x) Un{y), (7) 
n 

where the c„ are arbitrary constants, and the m„ are harmonic oscillator eigenfunctions 
corresponding to energy i5„ = {n+^)h. Obviously, this state is time independent. Nothing 
moves. If we try to get a semi-classical solution by using large values of n, we find that 
the amplitude of the wave function is large in the vicinity of the four corners of a square, 
x, ?/ ^ ±y/2En, and it is of course fixed in time. 

On the other hand, there is no such difficulty with the Heisenberg equations of motion, 
for example, dx/dr = [x,HT]/ih. These are formally identical to the classical oscillator 
equations of motion, and they lead to a nontrivial motion of the Heisenberg operators. 
We thus see that the Ehrenfest theorem [17] and, more generally, the correspondence 
principle, are not valid for such a dynamical system. 

In order to find a quantum counterpart to the dynamical system that is represented, 
in classical physics, by the Hamiltonian (|^), we must proceed more carefully. One of the 
harmonic oscillators, for example y, will serve us as a clock, and then the other one can 
be quantized in the usual way. We thus perform, still at the classical stage, a canonical 
transformation from y and Py to new canonical variables, 

= t^n-\py/y), (8) 

Po = -{pl + y')/2- (9) 

There is no corresponding unitary transformation in quantum mechanics (since the spec- 
trum is not invariant), but in classical mechanics, such a canonical transformation is 
perfectly possible. 

It is easily seen that [Q^, Hr]p^ = 1, so that 

dQ^dr = 1. (10) 
We can thus write 

H, = Po + H = 0, (11) 

where H = ^ (p^ + x^) is the ordinary Hamiltonian of the x-oscillator. Its equations of 
motion are dx/dr = and dpx/dr = —x. Thanks to Eq. (0), we can also write them as 

dx/dQ^ = px, (12) 

dpx/dQ^ = -X. (13) 



Finally, if we replace Px by —ihd/dx, and Pq by —iTid/dQ^, as usual, Eq. ( |Tl| ) becomes 
the standard Schrodinger equation for a harmonic oscillator. 
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However, at this point, we must be careful: the wave function ip{x, Q'^) should be 
normalized according to 

J \ij{x,Q')\'dx = l, (14) 

without any further integration / • • • dQ'^. This follows from our decision of keeping the 
clock time classical, so that it can play the ordinary role of time in Schrodinger's 
equation. In this way, wc have obtained a simple, consistent formalism, with the correct 
classical limit. Obviously, there are many other possible consistent formalisms, that are 
not equivalent to each other, and yet give the same classical limit. Quantization is possi- 
ble, but it is not a unique process. 



3 Definition of a dynamical time 

It will be now be shown how a similar quantization process can be performed for any 
classical dynamical system with a constrained Hamiltonian, 

Hr{q,p) = 0. (15) 

Here, q and p, without indices, mean {q^ . . . g"} and {pi . . -Pn}, respectively, and r is an 
arbitrary, convenient time parameter, in terms of which the problem has been formulated. 
In the case of Einstein's gravitational field equations, there is an infinite number of dy- 
namical variables and of constraints (there are twelve canonical variables, Qmn and tt"*", 
and four constraints per space point). In the present article, however, my main interest is 
in solving the "problem of time" for minisuperspaces. I shall therefore assume that there 
is only a finite number, n, of degrees of freedom, and a single Hamiltonian constraint, 
Eq. (15). 

Following the method illustrated in the preceding section, let us seek a canonical trans- 
formation from {g^} and {pk} to new canonical variables, {Q^} and {Pij}, with /i — 
0, . . . ,n — 1, such that 

dg7dr= [Q°,i7,]p3 = l. (16) 

It follows from Eq. (16) that 

Hr^Po + H{Q'' . . . Q'^-\ Pi . . . P„_i). (17) 

The latter equation defines an effective Hamiltonian H. Note that H does not depend on 
Pq so that, after we replace Pq by —ihd/dQ^, the classical equation Hr = becomes a 
Schrodinger equation for the new time and the {n—1) dynamical variables, . . . Q^~^ . 

The first step thus is to find a suitable clock time Q'^{q,p). This can easily be done, 
as least in a restricted domain of phase space, as shown in Fig. 1. We start with an 
arbitray (2n — l)-dimensional hypersurface /C, oriented in such a way that the fiow lines 
dq/dr = dHr/dp and dp/dr = —dHr/dq are nowhere tangent to /C. That is, all the 
fiow lines lie on the same side of /C, for r positive and short enough. Then, at least 
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(3^ = 



FIG. 1. Orbits in phase space, starting on the hypersurface /C. 
Each one of these orbits defines an internal clock-time, Q^{q,p). 

for some finite time, these fiow fines will not intersect — as long as they do not reach a 
critical point — and they will not reenter /C from the other side (however, if the motion is 
bounded, for example if it is periodic, reentry must obviously happen after enough time 
has elapsed). Anyway, for a finite time, each flow line that originates from /C ascribes a 
unique set of q and p to each value of r and, conversely, in a finite domain of phase space, 
there is a unique r for each set of {q,p}, say r = f{q,p). 

There still is here a formidable technical difficulty, because generic dynamical problems 
are not integrable: the number of constants of motion is usually less than the number of 
degrees of freedom, and the function f{q,p) defined above cannot be obtained in closed 
form and does not exist globally. For that reason, some authors [18] take the liberty of 
"fixing the gauge" by an arbitrary choice of the function f{q,p), leading to a form which 
is convenient for further work. It is not clear to me why this is permitted. This is also not 
necessary, because there do exist approximation methods for performing a sequence of 
canonical transformations which reduce the Hamiltonian to a normal form [19, 20]. These 
give approximate constants of motion, which are represented in phase space by "vague 
tori" and are useful for describing the dynamics over extended time periods. These tori 
remnants [21, 22] become important in quantum theory because, if their missing parts are 
small compared to 27rh, the quantum system behaves as if it were regular, with ordinary 
selection rules. 

Let us now return to the "problem of time." What we need is a canonical transforma- 
tion such that — f{q,p) is a prescribed function. As explained above, f{q,p) is chosen 
in such a way that 

[f{q,p),H,]p^ = l, (18) 
and therefore dQ'^/dT — 1. At this point, it is natural to ask what would happen if we 
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had chosen another initial hypersurface, say /C', leading us to a different time function, 
f'{q,p), say. We would then have 



[fiq,p)~fiq,p),H^U = 0, (19) 

so that (/ — /') has to be a constant of the motion. Either it is a function of H.^ or, if 
there are other, nontrivial constants of the motion, (/ — /') can be a function of them. In 
particular, if Hj^ = is the only constant of motion, (/ — /') can only be a mere number. 
Anyway, it does not matter for the sequel whether whether f{q,p) is uniquely defined, up 
to a numerical constant, or can be modified by adding a nontrivial constant of the motion 
(thus effectively giving a different version of the theory). 

Our task thus is is to find explicitly a canonical transformation that leads to the 
decomposition ([T7|). This can be done, in principle, by the solution of a first order partial 
differential equation of the same type as the Hamilton- Jacobi equation. It is easiest to 
use a generating function [23] of type Fi, that we shall write as S{q, Q). We have 

Pk = dS/dq\ (20) 

= -dS/dQr (21) 

Since S is time-independent (there is no explicit appearance of r in S*), the new Hamil- 
tonian is numerically equal to the old one, as in Eq. (|1^. 
To obtain S explicitly, we have to solve 

Q0.;(,.^). ,22) 

The various Q^, with > 0, are unspecified integration constants in the solution of 
Eq. (22). As in the Hamilton- Jacobi case, there is no guarantee that (22) has well behaved 
global solutions. However, it is always possible to achieve arbitrarily close approximations 
in a finite domain. An example is given in the next section. 

Once we have obtained S{q,Q), we get P^{q,Q) from Eq. (pT]). We can then invert 
these equations, in principle, and find q{Q,P). Likewise, Eq. (|20|) gives us p = p{q,Q), 
and since q{Q,P) is already known, this gives p = p{Q,P). All these results are then 
substituted in so as to obtain the explicit form of Eq. ([T7|) . Finally, that equation can 
be quantized in the usual way, replacing Q° by a new variable, t (recall that dQ^ /dr = 1), 
and Pq by —ihd/dt. Note, however, that the new parameter t is not a function of the 
space-time coordinates: it is a function, f{q,p), of the phase-space coordinates. This is a 
meaningful dynamical time, not a meaningless (gauge dependent) coordinate-time. 



4 Quantization of a minisuperspace 

Let us finally return to general relativity. As the simplest example, consider a spatially 
flat Friedmann-Lemaitre universe [24], with metric 
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The matter source is a massless scalar field 0, for which the energy density and pressure 
are 



P = (24) 

where natural units have been used: c = SttG = 1. 

The Einstein field equations, for the above metric and sources, become ordinary differ- 
ential equations for the three variables N{t), a{t), and In order to obtain a quantum 
version of this theory, the above differential equations must be obtainable as the Euler- 
Lagrange equations resulting from a Lagrangian [15]. It is easily found that a suitable 
Lagrangian, giving the correct equations, is 

L= {^a^^ -3aa'^)/N. (25) 
Note that does not appear in L, so that 

Pn = dL/dN = 0, (26) 
and therefore 

Pn = dL/dN = -L/N^ = 0. (27) 

The fact that L = is an initial value constraint imposed on the variables a, a, and 0. 
Likewise, we have 

p^ = dL/d(j) = 4)/N. (28) 
This is a constant of the motion, because dL/d(j) — 0. Finally, 

Pa = dL/da = -6 a a/N, (29) 
and 

Pa = dL/da = 3 (i 0^ - d^) /N. (30) 

As in Sect. 2, it is convenient to introduce an auxiliary time r by means of dr = Ndt. 
We then have a new Lagrangian, given by L^dr = Ldt, so that the action remains 
invariant. Furthermore, it is convenient to introduce, instead of the radial scale variable 
a, a new variable, v{t) = a^{t), which scales the volume element. We then have 

from which we obtain 

P0 = 1^ (32) 
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The corresponding Hamiltonian is 



Note that both Lr and Hr vanish weakly, as a consequence of (pTl) . The non-essential 
dynamical variable N{t) has thus been eliminated, but it has left a remnant, which is 
the initial value constraint, Hj^ = 0. The equations of motion resulting from the new 
Hamiltonian are: = const., Eq. (P5|), and 



dpy _ dHr _ Pd,^ Spy"^ 

dr ~ dv ~ 2v^^ 4 ■ ^^^^ 

Our task now is to find a dynamical time function, = f{v,Pv,P(f)), such that 
dQ^ /dr = 1. (Obviously, / is not a function of the cyclic variable 0, since the latter 
does not appear explicitly in the equations of motion.) In other words, we want a func- 
tion f{v,py,p^) that satisfies 

[f{v,py,p^),Hr]p^= 1. (36) 

For this, we have to solve the equations of motion explicitly. 

Substitution of (|33|) into the right hand side of (|3^) gives, after some rearrangement, 

l{dv/dTf + lH,v = lp/. (37) 

This looks like the elementary energy equation for free fall of a particle of unit mass, 
height V, and total energy fp^^, in a gravity field g = f-f^r- The solution is 



v = -IH,t'±^Ip^t, (38) 

where the integration constant was set so that v = when r = (in other words, the /C 
hypersurface is given by f = 0). Since by definition f > 0, the ± sign in (38) has to be 
the same as the sign of of p<^r. Note that we are not allowed to set Hr = at this stage: 
consistency of the method that was proposed in the preceding section requires that the 
equations of motion be valid for the entire phase space, not only for the orbits with initial 
conditions that satisfy H^- = 0. 



It is possible to solve directly (38) for r, and then to substitute (p^ in the result. 
However, it is simpler to proceed as follows. From (38), we have 



dv 

dr 



■lHrT±Jlp^, (39) 



whence, thanks to Eq. 



dv I 

vpv = -l-^ = H^TT ^Ip^- (40) 
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Thus, (BB) becomes 



V 



■lT{HrTT \JIp4>) = -zr{vp^T\/'ip<t,), (41) 



and therefore 

r = fiv,p^,,p^) = — — j^. (42) 



It is easy to verify directly that Eq. (|36D indeed holds. 

The next step is to find explicitly the transformation from the original canonical coor- 
dinates to the new ones, that include Q° and Pq. We have, from Eqs. ( ^2]) and (|4^), 



qo = - (43) 

-lv{dS/dv)±^{dS/d^y 

where S = S{v,(l),Q^,Q^). An obvious way for obtaining a solution is to separate vari- 
ables, namely, 

S = <f)Q^ + S'{v,Q'',Q^), (44) 
so that 

Q' = P^. (45) 
Rearranging Eq. (|^), we obtain 

1 _ ,ds' r,Q' 



whose solution is 



± VI — . (46) 



(47) 



We thus have 



dS' _^ ( 1 r^pA 

in agreement with (|4l|) . 
Note that 

P, = -dS'/dQ' = -\v/m\ (49) 
whence 

v = -lP,{Qy. (50) 
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When these equations for v and p„ are substituted into 



we obtain 



The reduced Hamiltonian if, defined by Eq. ([T7|), thus is 



H 



(51) 



(52) 



Recall that the ± sign in H is the same as the sign of p^/Q^ . Note that if Q° is considered 
as equivalent to the time r, the number of degrees of freedom has been reduced by 2: the 
variable N disappeared in the transformation from t to r, and the v and variables have 
been absorbed in the dynamical definition of a "clock-time" . 

Here, we must be careful and avoid expressing , in Eq. (52), by means of the right 
hand side of (E2|). This would give 



H 



[wrong] 



(53) 



Such a way of writing the Hamiltonian is not correct: it would give the true equations of 
motion for v and only if the initial conditions are set in such a way that Hj^ = in 
Eq. (|34D , namely 

= I {vp^f. (54) 
Indeed, we have from (p3D 

dv/dr = [v, H]p^ 
and 

dp^/dr = [p^, H]p^ = p^^/v^ 

and these agree with Eqs. (^) and ( ^7|) only if (|54D is satisfied. Therefore H in Eq. (|53| ) 
is not a valid, unconstrained Hamiltonian for this problem. Only H given by (|52| ) is 
acceptable. (More generally, the reader may easily verify that Pq = Hr — H has vanishing 
Poisson brackets with all the canonical variables only on the hypersurface Hr = 0.) 

We thus remain with the reduced Hamiltonian (|5^), and we may now replace in it 
by T. The only nontrivial equation of motion is 



[wrong) 



[wrong) 



(55) 



(56) 



dr 



PB 



± 



(57) 



whence (p — 
Schrodinger equation 



± ^/|lnr. Quantization is trivial: the wave function ijj{p^,T) satisfies a 



(5^ 
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so that 



^ = F{p^) exp (^^^ yfp^ In rj , (59) 

where F{p^) is an arbitrary function that takes care of normahzation. 

Obviously, only a superspace with a larger number of degrees of freedom can give an 
interesting theory. Unfortunately, "interesting" also means "nonintegrable" : the function 
f{q,p) is not in general well behaved (it is not "isolating") and approximation methods 
must be used [19-22]. 

Finally, the question must be raised whether the notion of a minisuperspace is a valid 
approximation for studying quantum gravity [25] . The arbitrary imposition of symmetry 
constraints on the gravitational field freezes almost all its dynamical degrees of freedom 
in a way that appears to be incompatible with the existence of quantum fluctuations. 
A similar dilemma arises in elementary classical mechanics, when we impose mundane 
mechanical constraints, such as restricting the motion of a mass to a two-dimensional 
surface. Classically, such a system is well defined. However, its quantization is not unique 
and it essentially depends on the nature of the constraining forces [26]. I hope to return 
to this problem in a future publication. 
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